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Abstract. We present follow-up observations of the mid-Infrared dark clouds selected from the ISOGAL inner 
Galaxy sample. On-the-fly maps of 13 CO, C ls O and the 1.2 mm continuum emission were conducted at the IRAM 
30-m telescope, showing spectacular correlation with the mid-IR absorption. The dark clouds are distributed as far 
as the prominent molecular ring at a distance of 3 to 7 kpc from the Sun. The clouds exhibit shapes ranging from 
globules to thin filaments down to < 1 pc in size. The on-the-fly images obtained in 13 CO and C ls O confirmed 
that the cores are dense, compact molecular emitters, significantly more massive than local dark clouds (more than 
1000 Mq) and lie within low activity Giant Molecular Clouds (GMCs). Ratios of the emission in the J = (2 — 1) 
and (1 — 0) transitions of 13 CO and C ls O show a remarkable uniformity within each cloud, with a significant 
portion of the sample represented well by a ratio of 0.67±0.12. Preliminary analysis of temperature and density 
measurements reveals that most of the cores have densities above 10 5 cm -3 and temperatures between 8 and 
25 K, these latter clouds being associated with young embedded stars. Despite the high extinction inferred from 
mid-IR (A v > 50, Hennebelle et al. 2001), the molecular lines are surprisingly weak, indicating likely depletion 
onto cold grains. 
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1. Introduction 

Galactic dark clouds observed in the mid-Infrared were 
first surveyed by ISO (Perault et al. 1996), and turn out to 
\ be the highly condensed parts of Giant Molecular Clouds 
(CMC's) kiloparsecs away from the Sun. The ISOGAL 
survey imaged ~ 10% of the inner Galactic ridge, mostly 
towards relatively quiescent areas, in broad filters around 
7 and 15 /xm (Omont et al. 1999). A systematic analysis 
of the ISOGAL plates (Hennebelle et al. 2001, hereafter 
Paper I) allowed extraction of a catalogue of about 450 
objects, most of them located in the inner Galaxy. The 
features are associated with GMCs lying between 3 and 
7 kpc from us. In Paper I, we derived opacities at 15 /xm 
in the range 1 to 4 for a few selected objects, leading to 
column densities of the order of 10 23 cm~ 2 . 

Similar findings have been reported by Egan et al. 
(1998) from the MSX (Midcourse Space Experiment) sur- 
vey of the whole Galactic plane. They counted about 2000 
IR dark clouds (IRDC's) located in the same distance 
range as the ISOGAL clouds and estimated extinctions 
at 8 fxm in excess of 2. Detections of few clouds at mil- 
limetre wavelengths (Carey et al. 1998) confirmed that 
these objects are dense (nn 2 > 10 5 cm~ 3 ) and cold (T 
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< 20 K). The authors conclude that these clouds are pre- 
protostellar cores where no sign of star-formation has been 
detected so far. This analysis was recently refined using 
continuum detection at 850 and 450 /im (Carey et al. 
2000). The relatively high masses inferred from the sub- 
millimetre observations suggest a significant potential for 
stellar formation in these cores. 

In the present paper we analyse spectroscopic and con- 
tinuum follow-up observations of the ISOGAL dark fea- 
tures, conducted at the IRAM 30-m telescope. A sample of 
13 objects has been mapped in several molecular tracers. 
A few characteristic physical parameters are inferred. The 
spatial resolution of these observations (11" at 1.3 mm) 
opens the way to the analysis of the gas associated with 
the dense dust revealed by the mid-IR absorption at an in- 
termediate spatial scale. Our main purpose is to compare 
the properties of our objects to those of the well-known 
local dark clouds. A subsequent study will provide an ex- 
tensive analysis of temperature and density measurements 
and better assess the physico-chemical processes at work 
in the clouds (Teyssier et al. 2001b). 

Sect. H presents the observational conditions and strat- 
egy. In Sect. [| we describe the correlation between the 
millimetre emission and the IR absorption data. The na- 
ture of the different objects is assessed and a detailed spa- 
tial and spectral analysis is introduced. Sect. || gives the 
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Table 1. Line parameters Table 2. Source parameters 
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(GHz) 
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(J2000.0) 


(J2000.0) 


(kpc) 


i ,j > ^ ^ 
"CO 


1^0 


110.201353 


22.5 


150-250 


DF+04.36-0.06 


17:55:53.07 


-25:13:18.7 


3.5 




2^1 
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11.2 


400-900 


DF+OQ 86-0 04' A ' 


1 8-07-37 99 
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c 18 o 


1 -> 
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DF+15.05+0.09 (A) 
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2 -> 1 
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11.3 
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DF+18.56-0.15 
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HC 3 N 
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81.881468 


30.2 
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DF+18.79-0.03 
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3.6 




10 -> 9 


90.978993 


27.2 


130-250 


DF+25.90-0.17 


18:39:10.13 
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5.5 




11 -■> 10 


100.076389 


24.7 


130-250 


DF+30.23-0.20 (s) 


18:47:13.16 
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6.7 
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5 -* 4 


85.457300 


28.9 
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DF+30.31-0.28 


18:47:39.03 


-02:27:39.8 


6.3 




6o — ► 5o 
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the rest frequency of the line, HPBW is the ideal half 
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and A the line wavelength, T sys 


is the SSB sysi 
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DF+36.95+0.22 


18:57:59.51 


+03:40:33.3 


5.0 












DF+51.47+0.00 (B) 


19:26:12.74 


+16:26:12.6 


5.3<°) 



physical properties of the clouds and their relation to the 
mid-IR opacities estimated in Paper I. Sect. [| summarises 
our conclusions. 



2. Observations 

Observations of 13 IR dark clouds were conducted at 
the IRAM 30-m telescope between July 1998 and August 
1999. Two summer runs were dedicated to spectroscopic 
observations, while continuum mapping at 1.2 mm was 
performed during winter time. Line and source parame- 
ters are listed in Tables [j] and ||. 



2.1. Spectroscopic observations 

Spectroscopic observations were performed during sum- 
mertime under good to average weather conditions, in 
single-side band (SSB) mode, with receiver temperatures 
of about 90 K at 3 mm and 180 K at 1.3 mm. The data 
were calibrated to the T£ scale using the chopper wheel 
method (Penzias & Burrus, 1973). Discussion on the final 
brightness temperature scale Xj^ b adopted for our data is 
presented in appendix |A.l| . Using comparison with obser- 
vations of line calibrators (Mauersberger et al. 1989), we 
believe our absolute calibration to be better than 20%. 
The pointing accuracy was checked by repeated contin- 
uum scans across planets and strong quasars. The spec- 
trometer was an autocorrelator set to a resolution of 80 
kHz, yielding velocity channels of 0.24 km s _1 (0.1 km 
s _1 ) at 3 mm (1.3 mm). Additional details on the tech- 
nique used to correct for baseline spurious effects in the 



autocorrelator are given in appendix A. 2 



The maps were obtained with the Spectral Line On- 
the-Fly technique (Ungerechts et al. 1999) with a scanning- 
speed of l"/second (samples every 2 seconds) and a cross- 
scan sampling of 6". We used two classes of reference posi- 
tions. The targets were mapped using a close-by reference 
position that was then compared to a further one chosen 



^ ^ Observed in April'99 (continuum) 

(s ' Observed both in April'99 and August'99 

More than one value may be given in case several lines are detected 
on the line of sight 

' r) ' Distance of the tangent point is assumed 



from the Massachusetts-Stony Brook Galactic Plane CO 
Survey (Sanders et al. 1986) and estimated as sufficiently 
free of emission. Finally, whenever possible, map coverages 
in the orthogonal direction were achieved allowing use of 
the PLAIT algorithm developed by Emerson et al. (1988). 

Maps of 3' x 3' or 4' x 4' were obtained for 7 of the fields 
and narrow bands of 0.5' x 3' for the other ones. HCO + and 
HCN J = 1 — lines were measured in parallel with 13 CO 
and C 18 (HCO + and HCN data are not discussed here). 
Examples of integrated maps are given in Fig. [j] and ^. 
Noise levels (rms) reached are around 0.26 K (0.53 K) at 3 
mm (1.3 mm) per regridded pixel and spectral resolution 
elements of 80 kHz. Three of the targets were observed 
in a second run in HC3N and CH3CCH (see transitions 
111 Table |l|). HC 3 N maps of 1.4' x 1.4' were obtained to 
further allow for accurate comparison of lines at different 
frequencies by restoring all data to a common synthetic 
beam. Transitions of CH3CCH were observed in position 
switching, yielding a 20 mK rms level per spectral resolu- 
tion element. 

2.2. Continuum observations 

Five clouds were mapped using the MPIfR 37-channel 
bolometer array centred at ~1.2 mm (Kreysa 1992, 
Kramer et al. 1998). We used the on-the-fly mode for con- 
tinuum observations, consisting of dual-beam rasters with 
scanning velocities of 6"/sec or 8"/sec. The spatial sam- 
pling interval in elevation was larger than generally used at 
the 30-m telescope. With a vertical offset of 24" (instead 
of 4-5"), individual detector maps are imcfer-sampled but 
the co-added final map is fully sampled (Teyssier & Sievers 
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Fig. 1. Maps of CO isotopomers integrated emission for DF+09.86-0.04. Contours are 1 to 4 Kkms' 1 by steps of 0.5 
for C 18 transitions and 4 to 12 Kkms^ 1 by steps of 1 for 13 CO maps (T^ scale). Velocity interval is 16-20 kms -1 . 
The white box indicates the core position while the black one frames the envelope area used in the analysis. 



1999). In a relatively short time (about 20 minutes) we 
could cover areas of 7'x6' with homogeneous observing 
conditions and a small field curvature when re-projecting 
in the equatorial frame. All maps were repeated several 
times at different hour angles and with different wobbler 
throws (46", 56" and 78") in order to allow the restora- 
tion of all spatial frequencies (see Emerson et al. 1995, 
or Pierce-Price et al. 2001). This additional information, 
however, is not used by the IRAM data reduction software 
(NIC, Broguiere et al. 1999) that we have used to derive 
the images presented here. A straightforward restoration 
algorithm (Emerson, Klein & Haslam 1979) is applied in- 
stead, followed by zero-order baseline substraction and a 
systematic skynoise removal. Fig. |^ gives an example of 
two of these maps. 

The beam size at 1.2 mm was measured to be ~ 11" us- 
ing Uranus. The calibration was achieved through regular 
on-thc-fly and on-off observations on planets (Uranus and 



Mars). We estimated both relative and absolute calibra- 
tion to be within ~ 10%. The zenith atmospheric optical 
depth was found to be between ~ 0.1 and ~ 0.3, according 
to regular antenna tipping. 



3. Observational properties 

3.1. Identification of the clouds and spatial distribution 

The infrared absorption features coincide remarkably with 
the C ls O emission ridges (left and middle panels of Fig. [|), 
as well as with the continuum emission, although some- 
times to a smaller extent (right hand panels of Fig. 
In some of the clouds local peaks in the millimetre maps 
coincide with the mid-IR emission of embedded stars or 
stellar clusters seen by ISO. The identification of OH and 
CH3OH masers (Caswell et al. 1995, see Fig. |3|) associ- 
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13 
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13 
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1 



30.235 30.23 30.225 30.22 30.215 
1 



Fig. 2. Same as Pig. |] for the central area of DF+30. 23-0.20. Contours arc 4 to 12 Kkms^ 1 by steps of 1 for C 18 
transitions, 16 to 35 Kkms" 1 by steps of 2 for 13 CO(J = 2-1) and 16 to 45 Kkms" 1 by steps of 3 for 13 CO(J = 1-0) 
(T^ scale). Velocity interval is 103-107 km s _1 . The black box indicates the filament position while the white one frames 
the star edge area. 



ated with these bright sources give additional evidence 
that stars may have formed recently. 

Fig. | and | illustrate such cases. For DF+09.86-0.04, 
the embedded stars appear as two small black squares 
indicated by arrows in the infrared data (marked with 
white stars in the continuum map). For DF+30. 23-0. 20, 
the stars themselves have not been mapped in CO, but 
are clearly seen as strong spots in continuum (white stars 
symbols). The dark filament turns to a bright filament 
in the mid-IR, while its millimetre emission is strongly 
reinforced around newly formed stars (Fig. |^). 

Some of the clouds present several emission line com- 
ponents along the line of sight (e.g. DF+30. 36+0. 11 and 
DF+31.03+0.27, which are the most distant clouds of our 
sample). The structures associated with different veloc- 
ity components present very different shapes and orienta- 
tions, which sometimes need to be combined to reproduce 
the shape of the dark cloud. For all clouds the velocity 
field exhibits significant structure, as illustrated for one 
component of the DF+31.03+0.27 cloud by the channel 
maps of Fig. [|. In this object, the profiles are complex 



and the gas structures seem to be connected over a wide 
velocity range. 

Large field maps obtained with the 4-m Nanten tele- 
scope in 13 CO(J =1 — 0) (Zagury et al., yet unpub- 
lished data) and shown on Fig. | for DF+09.86-0.04 and 
DF+30. 23-0. 20 indicate that the IR dark clouds are not 
isolated objects: they are embedded in quiescent CMC's, 
most of them belonging to the cold GMC population de- 
tected earlier in CO galactic surveys (e.g. Sanders et al. 
1986). This contradicts Egan et al.'s (1998) claim that the 
IRDCs are isolated objects. 

Using the galactic rotation model of Burton et al. 
(1991), we are able to derive kinematic distances of these 
objects from the Doppler shift of the line. Assuming that 
the kinematic distance ambiguity is solved by the absorp- 
tion bias in favour of the nearest position, we derive dis- 
tances between 2.8 and 7.4 kpc (see Table |^). For all ob- 
served clouds the inferred kinematic distances are consis- 
tent with objects lying within or in front of the molecular 
ring. 
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Fig. 3. Continuum maps at 1.2 mm. Left: DF+09.86-0.04. Contours are 10 to 60 mJy/11" beam by steps of 10. Right: 
DF+30.23-0.20 and surroundings. Contours are 10 to 100 mJy/11" beam by steps of 10, then 100 to 200 by steps 
of 50. CH 3 OH masers detected in previous surveys (Caswell et al. 1995) are indicated. White star symbols indicate 
embedded stars identified in the mid-IR plates. Positions are in galactic coordinates. 



3.2. Line intensities and velocity structure 

Figs. 7 and |§| give examples of the line profiles observed 
in three of the clouds for so-called core (C 18 emission 
peaks) and envelope (core periphery) areas. The lines, 
although often centrally peaked, rarely exhibit symmet- 
ric Gaussian profiles. In the 2 cases where embedded 
stars are associated with the clouds, the profiles appear 
as self-reversed in 13 CO(J = 2 — 1) and flat-topped in 
13 CO(J =1-0). DF+30.23-0.20 illustrates one of these 
cases (upper panel of the second column of Fig. 7). The 
line broadening due to a higher opacity in the star vicinity 
is quite significant. 

The absolute line intensities are smaller than what 
would be expected from such massive dense clouds if 
extrapolated from values observed in local dark clouds 
(Table ||). This effect is even stronger for weaker lines 
like HC3N or CH 3 CCH where the lines are a factor of 
5 weaker than in e.g. TMC-1 (Pratap et al. 1997). Part 
of this weakness may be attributed to a very low kinetic 
temperature or a small filling factor. But chemical effects 
or, most likely, depletion onto grains should be also con- 
sidered. On the other hand, the lines are broader than in 
local dark clouds: this is not surprising given the large size 
and mass of the objects. 

Line intensities ratios are compared in the scatter plots 
of Fig. § and|l0[ The 13 CO( J = 1 - 0) to 13 CO( J = 2-1) 
(noted R 13 ) and C 18 0(J = 1 - 0) to C 18 0(J = 2-1) 
(noted Ris) ratios are found to be constant on average 
for a given source, independent of the choice of pixels in 
the line centre or wings, in the cloud core or envelope. A 
similar result was found by e.g. Falgarone et al. (1998) 
for dense cores close to the Sun (d<150 pc). A signifi- 



Table 4. J = 2- ltoJ=l-0 ratios for C 18 and 13 CO 
estimated over the entire clouds. 



Name 


Ris 


Rls (A) 


Ris 


Rl 3 (A) 


DF+04.36-0.06 


0.66 


0.84 


0.59 


0.60 


DF+09.86-0.04 


0.96 


1.26 


0.93 


0.92 


DF+15.05+0.09 


0.62 


0.82 


0.72 


0.78 


DF+18.56-0.15 






0.59 


0.61 


DF+18.79-0.03 


1.02 


1.18 


0.66 


0.70 


DF+25.90-0.17 




0.70 


0.53 


0.60 


DF+30.23-0.20 (c) 


1.58 


1.63 


1.02 


1.01 


DF+30.36+0.11 


1.07 


1.17 


0.74 


0.87 


" (2 nd component) 


0.90 


1.24 


0.94 


0.79 


DF+31.03+0.27 


0.63 


0.85 


0.61 


0.69 


" (2 nd component) 


0.84 


1.07 


0.75 


0.82 


DF+51.47+0.00 


0.70 


0.78 


0.60 


0.63 


Mean 


0.90 


1.05 


0.72 


0.75 


(1-ct deviations) 


±0.28 


±0.27 


±0.15 


±0.13 



' ' Ratio of integrated intensities 

Dispersion is too large 
(c) Measured between 103 and 110 km s _1 . 



cant departure from this general behaviour is however ob- 
served in the close vicinity of the young stellar objects of 
DF+30.23-0.20: in the line wings the R13 ratio is higher 
than the average value obtained for the filament (right- 
most data points of Fig. [l(]). The excitation is certainly 
higher and the (J = 2 — 1) transition is optically thicker 
than the (J = 1 — 0) one, which results in a broader line 
and eventually in some self- absorption (Fig. 7). 

Table f| compiles the best-fit results of these ratios. 
The ratios between integrated intensities are also given 
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Fig. 4. Upper panels: Left: Contour map of the C 18 0(J =1 — 0) transition integrated emission for DF+09.86-0.04 
superimposed over the corresponding ISOGAL field in the LW3 (15 fan) filter. ISOGAL data have been reversed so 
that absorption features appear with positive values (see text). Arrows indicate positions of embedded stars. Middle: 
same for the C ls O(J = 2 — 1) transition. In both cases, contours are 1 Kkms -1 to 4 Kkms -1 by steps of 0.5 (T^ 
scale), velocity interval is 16 to 20 kms' 1 . Right: Contour map of 1.2 mm continuum emission for DF+09.86-0.04 
superimposed over the corresponding ISOGAL field in the LW3 filter. Contours are 10 to 60 mJy/11" beam by steps 
of 10. Lower panels: same as upper panels, for DF+30.23-0.20. Contours are 4 to 12 Kkms -1 by steps of 1 for both 
C 18 maps, and 10 to 170 mJy/11" beam by steps of 10 for the continuum map. Coordinates are galactic longitude 
and latitude. 



for comparison with previous 12 CO studies. These results 
can be compared to values observed in other clouds of 
our Galaxy: in quiescent clouds of the solar neighbour- 
hood, Falgarone et al. (1998) found a mean of 0.65±0.15 
for 12 CO and 13 CO. The C 18 observations reported by 
Kramer et al. (1999) in IC5146 lead to a mean value of 
0.83±0.22. Several studies have been conducted in the 
Galactic plane in 12 CO, showing a large-scale gradient 
of the integrated intensity ratios: 0.75 at ~ 4 kpc from 
the Galactic centre, to 0.5 at 8 kpc from the Galactic 
centre (e.g. Sakamoto et al. 1995, who derive a mean 
value of 0.66±0.01). A similar study conducted by Oka 
et al. (1998) in the Taurus and the Orion A Molecular 
Clouds, including star-forming sites, lead to mean values 
of 0.53±0.01 and 0.75±0.01 respectively. 



Our mean values of 0.72±0.16 ( 13 CO) and 0.90±0.29 
(C 18 0) are higher than the values found in these published 
studies. Even the starless sub-sample shows rather high 
values: 0.67±0.12 ( 13 CO) and 0.79±0.22 (C 18 0). Provided 
relative calibration uncertainties can be excluded, these 
values may be characteristic of the advanced condensation 
stage of IRDC's. 

The comparison of 13 CO and C 18 line intensities 
(called R 10 and R 2 i in Table [l] for J = 1 - and J = 2 - 1 
respectively) illustrate the differential abundances, excita- 
tion states and optical depths of the two CO isotopomers. 
At low intensities (yet above 5cr), these intensity ratios 
are consistent with the local interstellar [ 13 CO]/[C 18 0] 
ratio of about 7.3 (deduced from [ 16 0]/[ 18 0] = 560 and 
[ 12 C]/[ 13 C] = 77, Wilson & Rood 1994), but they depart 




from this value as intensities increase. The ratios drop 
below 2 close to the emission peaks. These ratios are sys- 
tematically larger in the envelope than in the core areas. 
The slow apparent saturation of the stronger line should 
be accounted for by detailed models. 

4. Cloud properties 

Preliminary results are derived from a simple analysis of 
the available data: spectroscopic intensities are interpreted 
with "LVG" calculations, which should be seen as the low- 
est order approximation (1-zone) to radiative transfer cal- 
culations. 

4.1. Temperatures and densities from spectroscopic 
data 

We used the spectra obtained in the J — 5 — 4 and 6 — 5 
transitions of CH 3 CCH (Fig. ||) to derive the kinetic tem- 
peratures in a subset of lines of sight. The details of this 
analysis will be presented in a subsequent paper (Teyssier 



et al. 2001b). The temperatures inferred are in the range 
8-25 K. The highest temperatures are representative of 
clouds with embedded stars (e.g. DF+30.23-0.20). These 
are also the only areas where significant emission is de- 
tected in the K = 2 rotation level. 



We use these temperatures to derive densities and 
molecular column densities from LVG simulations of 
HC 3 N, 13 CO and C 18 0. The HC 3 N transitions are rela- 
tively close in energy so that the error bars are somewhat 
large. The densities derived are larger than 10 5 cm~ 3 in 
the densest parts (cores). Observations of lower frequency 
transitions are included for more accurate derivations in 
our subsequent paper (Teyssier et al. 2001b). The corre- 
sponding masses vary between 2 x 10 2 and 2 x 10 4 M Q . 
The temperatures and densities of our sample clouds are 
similar to the parameters reported by Carey et al. (1998) 
for the MSX IRDC's. 




Fig. 6. Large field maps of DF+09.86-0.04 (left) and DF+30.23-0.20 (right) in the 13 C0( J =1-0) transition. Beam 
HPBW is 2.5' and sampling interval is 1'. Velocities are integrated between 16 and 20 km s _1 and 104 and 107 
km s _1 respectively. The spatial connections with larger scale structures is clear. White boxes indicate the extent of 
the bolometer maps displayed in Fig. [| 



4.2. Column density comparison 

For several of the sources we can compare independent 
estimates of total column densities: from mid-IR opac- 
ities, from dust 1 mm emission, from densities and size 
estimates, or from the column densities of tracers. 

The first method consists in relating the mid-IR opac- 
ities to visible extinctions, hence total gas column density. 
The dust emissivity model of Draine & Lee (1984, here- 
after DL84) implies: 



A v ~ 65 x r 7/ir 



(1) 



The total column density along the line of sight can then 
is estimated using the (Nh + h 2 ,Av) relation of Bohlin et 
al. (1978). 

The second method is based on continuum measure- 
ments of the optically thin dust at 1.2 mm and uses the 
relation described in Motte et al. (1998): 



Si, 



(2) 



where S mm is the 1.2 mm flux density in mJy per 
11" beam, \i the mean molecular weight, mn the mass 
of atomic hydrogen, Bi.2(T) the Planck function at wave- 
length 1.2 mm and temperature T, and k±.2 the dust mass 
opacity at 1.2 mm. We adopted a value — 0.003 

cm 2 g _1 , consistent with the DL84 dust model and the 
(-/Vh+h 2 >Av) relation used in method (1). Such a value 
is believed to be suited to pre-stellar dense clumps and 
cores with a typical uncertainty of a factor of ~ 2 (see 
Motte et al. 1998 and ref. therein for a complete discus- 
sion). The dust temperature is assumed to b e equal to the 
gas temperature as estimated in section 4.1. 



In the third method we assume that the H2 density 
derived from beam-averaged spectra is representative of 
the density within a sphere of projected size in the sky 
equal to the beam diameter. This assumption seems vi- 
able, in view of the apparent size of the core regions, but 
might be very uncertain, if the gas distribution is lacunar. 
Estimates are given both from HC3N (method (3)) and 
from C 18 (method (3')). 

Finally (method (4)), we used the (N K2 ,LF(C 18 0)) re- 
lation derived by Cernicharo & Guelin (1987, hereafter 
CG87) in local dark clouds. Table || gathers our results. 

The column densities range from a few 10 22 to at most 
a few 10 23 cm -2 . The ratio of the estimates from methods 
(1) and (2) is within ± 15% of 0.65 for all 4 cases where 
both estimates are available. A higher value of the emis- 
sivity at 1.2 mm (1.6 times the DL84 value) is required 
for matching the 2 derivations. This agrees with obser- 
vations of cold condensations in the solar neighbourhood 
(e.g. Dupac et al. 2001). A proper account of the dust tem- 
perature distribution is in all cases requested for deriving 
more accurate conclusions. 

Estimates from method (4) are 2 to 9 times smaller 
than values traced by the dust, the better agreement cor- 
responding to the embedded star cluster of DF+30.23- 
0.20. Since most of the reported points have high A v , this 
trend suggests possible molecular depletion onto very cold 
grains (see next section). 

Since we considered the same core size for both cases 
(3) and (3'), the discrepancy between these two estimates 
is not surprising: it illustrates how the density derivations 
are limited by the different critical densities (e.g. Evans 
1980). On the other hand, the column densities derived 
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Fig. 7. Average spectra of CO isotopomers for DF+09.86-0.04, DF+30.23-0.20 and DF+31.03+0.27. First row of 
panels: 13 CO(J =1 — 0) (full lines) and 13 CO(J = 2 — 1) (dashed lines) transitions in the core area. Second row 
of panels: C 18 0(J =1—0) (full lines) and C 18 0(J = 2 — 1) (dashed lines) in the core (or star vicinity) area. Third 
and fourth rows of panels: same as 1st and 2nd rows for a different area of the same clouds. The spectra have been 
convolved to identical resolution (22"and 0.213 km s -1 channels), and averaged over 30" boxes. Temperatures are in 
the T^ h scale (see appendix A.l). 
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Table 3. CO line parameters for spectra averaged over 30" x 30"areas 
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In each case, we give the peak temperature (K, in T^ b scale, see appendix A.l) with its noise r.m.s. (lcr r1Tls ) in brackets and the line width 
(kms - , FWHM — 2.35<r). Temperatures correspond to data smoothed to the same spatial (22") and velocity (0.213 km s - ) resolution. 
Brackets give the uncertainty on ratios based only on the spectra r.m.s. noise. 
Extrapolated from self-reversed spectra. 

Overlap with star clusters identified on the ISOGAL plates 



from HC3N ((3)) are generally higher than those inferred 
from the dust ((1),(2)). This may indicate that the actual 
volume filling factor is smaller than 1, the value assumed 
for this estimate. As in the picture proposed by Lada, 
Evans & Falgarone (1997) in local dense cores, this result 
suggests that the densest regions of the clouds are further 
fragmented. 

4.3. Extinction 

Using Eq. [j], we studied the pixel-to-pixel correlation be- 
tween visual extinction and the C 18 0(J =1 — 0) inte- 



grated intensity for selected clouds. The error estimate on 
the mid-IR opacities is described in appendix B. Fig. [Tl| 
shows the resulting scatter diagrams. The scaling of ex- 
tinction data to A v depends on the assumed dust emissiv- 
ity model. We then denote the visual extinction scale aA v , 
where a = 1 for the DL84 model. The result of bivariate 
fits of our data are displayed on Fig. [ll]. The errors result- 
ing from the linear regression show that this best fit line 
is relatively well constrained, regardless of its adequacy, 
or not, to properly represent the data over the whole A v 
range, of the intrinsic data dispersion and of the absolute 
calibration uncertainty. 
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Fig. 8. Sample spectra of HC 3 N, CH 3 CCH(J = 5-4) and CH 3 CCH(J = 6-5) for DF+30.23-0.20 in the star vicinity. 
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Fig. 9. Scatter plots of CO iso- 
topomer pairs of temperatures T' mh 
for DF+09.86-0.04 (see text). The 
data have been smoothed to a 
22" spatial and 0.2 km s _1 spectral 
resolution. Velocities considered are 
between 16 and 20 kms -1 . Each 3-d 
pixel is represented by a dot. Best- 
fit results are plotted for the upper 
panels. Local ISM ratios are indi- 
cated by dashcd-lines in the bottom 
diagrams. 



We now compare our results to earlier studies of local 
dark clouds, based on star counts. In HCL2, CG87 inferred 
for A v ranging from 1.5 to 6: 

W(C 18 0) = 0.28 ± 0.05 • (A y - 1.5 ± 0.3) K km s" 1 (3) 

For A v < 10, Alves et al. (1999) found in L977: 

W{G 18 0) = 0.18 ± 0.01 ■ (A w - 1.67 ± 0.26) K km s" 1 (4) 

These relations are overplotted on the graphs of Fig. [ll]. 
An important difference with the earlier studies is that 



the extinction ranges probed barely overlap. The different 
location and nature of the sampled clouds observed with 
different spatial resolutions might lead to large discrepan- 
cies. The different data sets nevertheless seem compatible. 
Our scatter diagrams exhibit a similar break as reported 
by Alves et al. (1999, see also Frerking et al. 1982, Lada et 
al. 1994) above A v ~ 10-15 mag. In our case, the depar- 
ture from the linear fit occurs at higher A v in the range 
15-25 magnitudes, which is about the sensitivity limit of 
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Fig. 10. Same as Fig. g for 
DF+30.23-0.20. The velocity 
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Table 5. H2 Column densities averaged over a 20" beam, in units of 10 cm 
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The density could not be constrained by the LVG model 



the previous studies. It is interesting to note that the re- 
lations given in Eqs. |^ and |] are not inconsistent with 
our few sample points at A v below 15-20 magnitudes (e.g. 
DF+18.56-0.15, DF+25.90-0.17). 

As for previous properties, our clouds are not all iden- 
tical and their behaviour in the (VF(C 18 0), A v ) plane may 
characterise different conditions: 

1. DF+15.05+0.09 and DF+18.56-0.15: these clouds 
are identified as the most opaque objects of the sample 



and thus suffer most from the mismatch in the extinction 
ranges. Nevertheless the low extinction relations (Eqs. || 
and |J) smoothly connect to our data points and are con- 
sistent with a break at A v around 15, similar to Alves 
et al. (1999) findings. Their conjecture that such a break 
is due to molecular depletion onto grains is validated by 
the Kramer et al.'s (1999) study based on rare CO iso- 
topomers. The similar break found here reinforces our as- 
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Fig. 11. Relation between the C 18 0( J =1 — 0) integrated intensity and the visual extinction derived from mid-IR 
measurements in six of the clouds, a depends on the A v /Ai5^ m scaling; a — 1 for DL84. IR data wer e sm oothed to 
the C 18 beam size of 22". The error in the extinction measurements is calculated according to Eq. B.5 while the 
uncertainties on W(C 18 0) correspond to the rms noise of the spectra. The solid line represents the least-square fit of 
the data weighted in both coordinates. The dashed (respectively dashed-dot-dot-dot) lines give the fits obtained for 
nearby dark clouds by Cernicharo & Guelin (1987) (resp. Alves et al. 1999) and recalled in Eq. || (resp. Eq|J). 



sumption of significant molecular depletion inside the IR 
dark clouds. 

2. DF+ 18.79-0.03 and DF+25.90-0.17: these clouds 
present much more points in the extinction range of Eqs. || 
and ^. Although the samples are relatively scattered, we 
find that the data are nicely fitted by the relation of CG 
and Alves et al. (1999) within 10-30%, which is still con- 
sistent given the uncertainty implied by the factor a intro- 
duced above. The break position is, again, consistent with 
previous results and the negative intercept in agreement 
with Alves et al. (1999) claim that shielding is probably 



not efficient enough to avoid molecular photodissociation 
from interstellar radiation at low visual extinctions. 

3. DF+9.86-0.04 and DF+30. 23-0.20: these clouds 
exhibit a clear lack of points of weak C 18 emission at 
low visual extinctions. As mentioned in previous sections, 
these objects are associated with several young star clus- 
ters whose radiation probably heatens the surrounding 
core. In DF+30.23-0.20, several OH and CH 3 OH masers 
(Caswell et al. 1995) have been identified at the edge of 
the filament (see Fig. |§|), probably tracing an even more 
evolved state of star-formation. 
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We conclude that our study extends earlier studies 
based on star counts in a consistent way. In particular 
it confirms that C 18 can be used as a rough column den- 
sity, and thus mass, tracer in dense dark clouds, provided 
corrections are made to account for progressive saturation. 

5. Conclusions 

In this paper we presented millimetre follow-up observa- 
tions of Galactic infrared dark clouds discovered in the 
galactic ISO survey at 7 and 15 /im. All objects are de- 
tected in 13 CO and C 18 and present remarkable spatial 
correlation with the mid-IR absorption data. We found 
that these clouds are not isolated but appear as the most 
condensed parts of quiescent GMC's 2 to 8 kpc away from 
us, some of them being associated with embedded stars 
or clusters observed by ISO. The clouds present a variety 
of shapes and line profiles indicative of supersonic flows. 
Velocity maps exhibit structures down to the parsec size 
connected to larger clouds along the line of sight. 

Special attention was paid to the intensity calibra- 
tion of our data and a simple approximation for extended 
sources seen in a main+error beam pattern was developed. 
Although not completely correct, it is believed to provide 
a better estimate of the main error beam integrated in- 
tensity than the widely used main-beam brightness tem- 
perature. 

In spite of a certain diversity, some common behaviours 
are observed in our cloud sample. We found that the 
(J = 2 - 1) to (J = 1 - 0) ratios of 13 CO and C 18 
were remarkably uniform within each object and that a 
significant amount of them match the values between 0.6 
and 0.8 observed in local dark clouds. Higher ratios are ob- 
served in the other objects, indicative of a broader range 
of conditions. 

A preliminary analysis of data obtained with temper- 
ature and density probes such as HC3N and CH 3 CCH 
shows that the clouds are mostly cold gas and dust con- 
densations of temperatures below 10 K and densities in 
excess of 10 5 cm~ 3 . In those conditions, depletion onto 
grains is expected to strongly affect the molecular emis- 
sion, but opacity effects can so far not be completely ruled 
out. Column densities close to 10 23 cm -2 are inferred from 
the dust absorption and emission. These clouds are in an 
advanced condensation stage, and some of them already 
contain young stars. Similar massive dark clouds are re- 
ported by Egan et al. (1998) and Carey et al. (1998) in 
their MSX survey of the Galactic plane: the parameters 
they derive are in good agreement with ours. 

The darkest condensations (Nn 2 > 10 23 cm~ 2 ) rep- 
resent ~ 1% of the inner Galaxy surveyed with ISO. 
Assuming that < 80% of the population escaped detection 
(distance bias, confusion,..), a wild extrapolation leads to 
at most 10% of the dust in the inner Galaxy trapped in 
these cold dense cores. This is a factor <~ 50 too low to ac- 
count for the hypothetical cold dust component advocated 
by Reach et al. (1995). 



Comparison of dust mid-IR absorption and 1.2 mm 
emission points towards emissivities at 1.2 mm higher (fac- 
tor <~ 1.6) than Draine & Lee (1984). On the other hand 
spectroscopic tracers indicate that the denser parts of the 
cold condensations are not volume filling, as earlier found 
in other situations (e.g. Lada et al. 1997). 

We studied the relation between the integrated C 18 
emission and the visual extinction and compared it to pre- 
vious results obtained for local clouds. Our study samples 
higher extinctions than these former studies but we find 
that the correlations consistently connect in the A v — 15- 
25 mag range and indicate that the C 18 emission still 
is a relatively good tracer of the cloud mass at visual ex- 
tinction above 15-20 magnitudes, provided corrections are 
made to account for radiative transfer saturation and/or 
depletion onto grains. 



Appendix A: Observation and reduction 
techniques 

A.l. Brightness temperature scaling 

When observing extended non-uniform emission, a signifi- 
cant fraction of the signal is peaked by diffraction pattern 
components larger than the main beam. As a consequence, 
the widely used main beam brightness temperature scale 
(T mb , also defined as by Kutner & Ulich, 1981) over- 
estimates the actual physical line intensity derived from 
the observations. Several methods have been proposed to 
correct from these effects (see Bensch et al. 2001a for a 
recent review) but most of them require additional ob- 
servations to be conducted on larger scale and generally 
smaller telescopes. 

In this work, we have used an approximation which 
consists of assimilating the source to a uniform disk-like 
emitter of brightness temperature T^. If the complete 
beam of the instrument is known, it is possible to cal- 
culate the ratio between T m b and the real 7b as a function 
of the source size. The 30-m beam is modelled according 
to Greve et al. (1998) as the superimposition of a main 
beam and three extended error beams: 

P = p i = f[ exp[-ln2 (26/9i) 2 } d6 (A.l) 

i=0,3 i=0,3 1 

where Pq is the main beam beam pattern and 

Wi = if Pidil (A.2) 

JJ2-K 

is the power fraction contained in the beam component 
Pi, and f2; its corresponding solid angle. Using the defini- 
tion by Kutner & Ulich (1981), the main beam brightness 
temperature can be written: 

jj s T h *Pdn 
Tmb " IL b p an (A ' 3) 
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For Gaussian beams and a disk-like uniform source of di- 
ameter 8 S and temperature Tb, Eq. [A. 3] finally is written: 



t —Th. 



i - A ex P [-in2 {e s /e,) 2 ] 



i=0,3 



with 



P dQ = Wi (Panis 1995) 



Similarly, one has: 

XX = /3o X mb , where /3 = 



X, 



(A.4) 



(A.5) 



(A.6) 



is the ratio of the main b eam to forward efficiency (e.g. 
Kramer et al. 1997). Table A.l lists the ratios for various 
disk diameters. As expected, effects at 3 mm remain neg- 
ligible if the source is smaller than the first error beam 
diameter. 

To estimate the accuracy of this crude approximation, 
we compared it to a more accurate correction used by 
Falgarone et al. (1998) and Bensch et al. (2001b) on their 
data. We found an agreement within 10% for the bright- 
ness temperatures at 1.3 mm, where the error beam con- 
tamination is rather significant, and within 7% at 3mm, 
where the error beam peak-up is much smaller. 

For our data, the corrections to T m b amount to 5-7 % 
at 3 mm and to values in the range 10-70 % at 1.3 mm. 
Emission at a larger distance from the observed position, 
which is not considered here, might increase the 3 mm 
correction a little. For most sources the adopted correc- 
tion at 1.3 mm corresponds to emission extending beyond 
the first error beam (~ 2' diameter). One of the greatest 
limitations of the method shows up for elongated sources 
but still the geometrical mean should represent on aver- 
age the fraction of the source surface picked up beyond 
the main beam. Despite this sometimes severe limitation, 
we conclude that our approximate correction reasonably 
accounts for error beam pick-up and that the residual er- 
ror on Tb is significantly smaller than the error resulting 
from using either T m b or T£. 

A. 2. Spurious effects in the spectral baseline 

Our spectral baselines are affected by standing waves (re- 
sulting from multiple reflexions along the optical paths) 
and by discontinuities in the autocorrelator (AC) baseline 
(gain mismatch between sub-bands, also known as plat- 
forming). The platforming can be efficiently corrected if 
a wider band backend (here a filter bank) is connected 
in parallel. This permits us to provide an accurate 0- level 
to each sub-band, independent of baseline ripples that are 
seen equivalently in both spectrometers. A sub-band read- 
justment based on 0-ordcr baselines separately estimated 
in each sub-band of the AC alone would indeed alter the 
wideband standing wave pattern and would thus bias any 
further correction of this additional effect. 



Table A.l. Error on brightness temperature for a disk- 
like uniform source of diameter # s 



B 


3 


mm 


1.3 


mm 


(") 


Tmb/Tb 


n/n 


Tmb/Tb 


n/n 


26.2 


0.67 


0.47 


1.00 


0.42 


30 


0.77 


0.54 


1.02 


0.43 


46 


1.00 


0.70 


1.06 


0.45 


60 


1.04 


0.73 


1.10 


0.46 


90 


1.05 


0.74 


1.20 


0.51 


120 


1.05 


0.74 


1.32 


0.56 


150 


1.06 


0.74 


1.44 


0.60 


180 


1.07 


0.75 


1.54 


0.65 


240 


1.08 


0.76 


1.68 


0.71 


600 


1.14 


0.80 


1.90 


0.80 


6000 


1.43 


1.00 


2.38 


1.00 



Appendix B: Error calculation for mid-IR opacities 

The error estimate on the opacities calculated in Paper I is 
based on the following formalism for the mid-IR emission: 



hot = h R e 



he 



(B.l) 



where h g and h g are respectively the background and 
foreground emission w.r.t. the cloud of opacity r. We as- 
sume that the main uncertainty on r comes from the in- 
trinsic fluctuations of hg and h g , of total rms amplitude 
e, which we estimate in unobscured neighbouring regions. 
Using the parameter (3 introduced in Paper I: 



hot 

the fluctuations of h g and h g are written: 



(B.2) 



(B.3) 



Assuming that these contributions are uncorrelated, this 
leads to: 



5t 2 



Vibg' 
^ hot ■ 



(— V 



l-(3e- T 
J 2 



(B.4) 
(B.5) 



This approach tends to favour the lowest extinction 
points of our samples but takes into account the fact that 
errors on nearby objects (/3 < 1) are smaller. 
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